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Abstract. For the well studied halo star HD 122563 and the four stars in the globular clus- 
ter NGC 6397, we determine NLTE abundances of iron using classical plane-parallel model 
atmospheres. Each star reveals a discrepancy in abundances between the Fe i lines arising 
from the ground state and the other Fe i lines, in qualitative agreement with the 3D-LTE 
line formation predictions, however, the magnitude of the observed effect is a factor of 2 
smaller compared with the predicted one. When ignoring the Fei low-excitation lines, the 
NLTE abundances from the two ionization stages, Fe i and Fe n, are consistent in each inves- 
tigated star. For the subgiants in NGC 6397, this is only true when using the cooler effective 
temperature scale of Alonso et al. (1999). We also present full 3D-LTE line formation calcu- 
lations for some selected iron lines in the solar and metal-poor 4480/2/-3 models and NLTE 
calculations with the corresponding spatial and temporal average (3D) models. The use of 
the (3D) models is justified only for particular Fei lines in particular physical conditions. 
Our NLTE calculations reproduce well the centre-to-limb variation of the solar Fe i 7780 A 
line, but they are unsuccessful for Fe i 615 1 A. The metal-poor (3D) model was found to be 
adequite for the strong Fei 5166 A (£ cxc = 0) line, but inadequite in all other investigated 
cases. 

Key words. Line: formation - Stars: abundances - Stars: atmospheres 



1. Introduction 

Most stellar-parameters and chemical- 
abundance studies of cool stars still employ 
classical plane-parallel (ID) model atmo- 
spheres. How large are the abundance errors 
caused by ignoring hydrodynamical phenom- 
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ena (3D effects) in the atmosphere of such 
stars? Collet et al. (2007) published (3D-1D) 
abundance corrections for some fictitious 
Fei and Fen lines in a small grid of model 
atmospheres representing the atmosphere of 
cool giants of various metallicities. These 
data were updated by Hayek et al. (2011) 
by including Rayleigh scattering in the line- 
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formation calculations. Significant changes 
were found only for the ultra-violet lines, with 
a wavelength of 3000 A. Therefore, for the 
lines lying in the visual spectral range, we cite 
hereafter the more extensive calculation of 
Collet et al. (2007). 

They showed that the 3D corrections are 
negative for the low-excitation (E exc < 2 eV) 
lines of Fe i independent of their strength in the 
[Fe/H] = -2 and -3 models. The magnitude 
of the correction depends strongly on excita- 
tion energy of the lower level. For example, in 
the r eff /log g/[Fe/H] = 5035/2.2/-2 model, the 
(3D- ID) correction is slightly negative for the 
£exc = 4 and 5 eV lines, but it goes down to 
-0.6 dex for the E exc = lines. The magnitude 
of correction increases towards lower metallic - 
ity. Classical abundance analysis of metal-poor 
(MP) stars should reveal an abundance differ- 
ence of 0.5 dex between different excitation 
Fe i lines if the theoretical predictions of Collet 
et al. (2007) are correct. 

The theory predicts positive corrections 
for the Fen lines. For example, (3D-1D) = 
+0.2 dex for the E &xc — 3 eV lines in the 
5035/2.2/-2 model. With 3D effects having 
opposite sign for Fei and Fen lines, rather 
large shifts to stellar surface gravities should 
result from spectroscopic analyses. This is 
worrisome given the overall decent agreement 
between spectroscopic analyses and stellar pa- 
rameters from stellar-structure models. 

This study carefully checks the iron lines 
in the few MP stars with well determined stel- 
lar parameters using classical ID model atmo- 
spheres and non-local thermodynamical equi- 
librium (NLTE) line formation. We are inter- 
ested in uncovering any discrepancies between 
Fe i and Fe n and between Fe i lines of different 
excitation energy and to evaluate their magni- 
tude if they exist. The obtained results are pre- 
sented in Sect. 3. 

At present, there are no tools to perform 
full 3D-NLTE line formation calculations for 
atoms with a complicated term structure like 
those of Fei and Fen. The full 3D line- 
formation calculations differ from their ID 
counterparts due to different mean atmospheric 
structures and the existence of atmospheric in- 
homogeneities. The effect of different mean at- 



mospheric structures on emergent fluxes and 
derived chemical abundances can be evaluated 
using the spatially and temporally averaged at- 
mospheric stratification denoted (3D). We in- 
vestigate here how closely the (3D) line forma- 
tion for iron resembles that of the correspond- 
ing full 3D model. Section 4 compares the ob- 
served solar centre-to-limb variation of some 
Fei lines with the predictions from the (3D) 
NLTE line-formation modelling. The LTE and 
NLTE calculations were also performed with 
the average 3D models representing the atmo- 
spheres of very metal-poor (VMP) cool gi- 
ants (see Sect. 5). Our conclusions are given 
in Sect. 6. 

2. Method of calculations 

The statistical equilibrium of Fe i-Fe n was cal- 
culated with a comprehensive model atom us- 
ing the method of Mashonkina et al. (2011) 
and the code DETAIL (Butler & Giddings 
1985), where the opacity package was up- 
dated. Inelastic collisions with H i atoms were 
taken into account employing the classical 
Drawinian rates scaled by a factor S H = 0.1. 
The line profile calculations were performed 
with the code SIU (Reetz 1991). 

For classical ID model atmospheres, we 
used the MARCS grid (Gustafsson et al. 2008) 
and the model 4600/1. 6/-2.5 computed with 
the code MAFAGS-OS (Grupp et al. 2009). 
The 3D model atmospheres were computed 
with the code C0 5 BOLD (Freytag et al. 2012). 
The average (3D) models were obtained by 
horizontally averaging each 3D snapshot over 
surfaces of equal (Rosseland) optical depth. 
For comparison, we also considered the ID 
hydrostatic mixing-length model atmospheres 
which employ the same micro-physics and ra- 
diative transfer scheme as C0 5 BOLD. They 
are referred to as 1D LH d models. 

In this study, we determined absolute stel- 
lar iron abundances based on the linelist pro- 
vided by Mashonkina et al. (2011). However, 
only the iron lines with accurate atomic data 
available were included in the analysis. For 
Fe i, we selected the lines with gf- values mea- 
sured by either the Oxford or Hannover group. 
The two different sets of gf- values were ap- 
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Fig. 1. The best fit (continuous curve) of the Fei 
5247 A line in HD 122563 (open circles) obtained 
from the NLTE calculations with the MAFAGS-OS 
model 4600/1. 60/-2. 5. The theoretical profile was 
convolved with a profile that combines instrumen- 
tal broadening with a Gaussian profile of 3 kms~' 
and broadening by macroturbulence with a radial- 
tangential profile of 4 kms~'. The differences be- 
tween observed and calculated spectra (O - C) mul- 
tiplied by a factor of 2 are presented in the upper 
part of figure (offset by +1.1). 
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Fig. 2. NLTE (Sh = 0.1) abundances derived from 
the Fei (open circles for £ cxc > 2 eV and rombs 
for E txc < 2 eV) and Fen (filled circles) lines in 
HD 122563 as a function of excitation energy of 
the lower level (top panel) and observed equivalent 
width (bottom panel). In each panel, the dashed line 
shows the mean iron abundance determined from 



plied for Fe n. One is of Melendez & Barbuy 
(2009, hereafter, MB09), and another is of 
Raassen & Uylings (1998, hereafter, RUm). 
The RUm ^/-values were corrected by a fac- 
tor +0.11 dex following the recommendation 
of Grevesse & Sauval (1999). 



3. 1D-NLTE analysis of iron lines in 
some selected metal-poor stars 

3.1. VMP cool giant HD 122563 

HD 122563 was selected as one of the best ob- 
served halo stars, with well determined stellar 
parameters. Recent measurements of its angu- 
lar diameter (Creevey et al. 2012) give T e s = 
4600+41 K. In our earlier paper (Mashonkina 
et al. 2011), we employed the same value, 
which was based on emprical colour calibra- 
tions, to calculate the gravity log g = 1 .60+0.07 
from the HIPPARCOS parallax. 

A high-quality spectrum of HD 122563 
{R - 80000 and S/N > 200) was taken from 
the ESO UVESPOP survey (Bagnulo et al. 
2003). The quality of line profile fits is illus- 
trated in Fig. 1 for Fe 1 5247 A. 



The NLTE abundances derived from indi- 
vidual lines in HD 122563, in total, from 29 
lines of Fei and 15 lines of Fen, are pre- 
sented in Fig. 2. We adopted the microtur- 
bulence velocity = 1.95 kms -1 as deter- 
mined by Mashonkina et al. (2011) from Fei 
and Fen lines. Interestingly, the use of dif- 
ferent sets of ^/-values leads to very sim- 
ilar mean abundances from the Fen lines, 
i.e., loge Fe ii(MB09) = 4.92 + 0.08 and 
logepeii(RUm) = 4.91 + 0.09. When calculat- 
ing the mean abundance from the Fe i lines, we 
included only those with E exc > 2 eV, which 
are predicted to be less affected by 3D effects. 
In NLTE (Sh = 0.1), we find consistent abun- 
dances between Fe i and Fe n, with the differ- 
ence (Fei(> 2 eV) - Fe ii)nlte = -0.06 dex, 
when using the MB09 ^/-values. In contrast, 
the difference between the corresponding LTE 
abundances amounts to -0.26 dex. We note 
that the Fe i lines arising from the ground term 
(.Eexc < 0.2 eV) reveal, on average, a 0.2 dex 
higher abundance compared to that of the re- 
maining Fei lines. 
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Collet et al. (2009) determined the LTE 
abundances from individual Fe i and Fe n lines 
in HD 122563 based on full 3D line-formation 
calculations. The stellar parameters employed 
were the same as in our study. For the Fei 
lines, Collet et al. (2009) found negative 3D 
abundance corrections, which depend strongly 
on the lower-level excitation potential. For ex- 
ample, (3D - ID) -0.3 for £ exc = 3 eV. In 
contrast, they obtained approximately 0.4 dex 
higher 3D than ID abundances from the Fen 
lines. Collet et al. (2009) argued that "such 
trends may be attributed to the neglected depar- 
tures from LTE in the spectral line-formation 
calculations". 

When applying the 3D corrections of 
Collet et al. (2009) to our NLTE abundances, 
we obtain a dramatic discrepancy between Fe i 
and Fen, with (Fei(> 2 eV) - Fe ii)nlte+3d = 
-0.75 dex. 

As for the Fe i lines arising from the ground 
state, 3D acts in right direction, but the pre- 
dicted magnitude -0.8 dex of the 3D correc- 
tion is markedly too large. 

3.2. Subgiants and red giants of 
NGC 6397 

Our analysis was extended to higher tem- 
peratures and higher gravities by inspecting 
the iron abundances of stars in the globu- 
lar cluster (GC) NGC 6397 on high-quality 
(R = 47 000, S/N 100) spectra taken with 
the FLAMES-UVES spectrograph at the Very 
Large Telescope (see Korn et al. 2007, for a 
detailed description of the observations). We 
selected the two stars N5281 and N8298 on 
the subgiant branch (SGB) and the two stars 
N11093 and N14592 on the red giant branch 
(RGB). Stellar parameters of these stars were 
discussed in detail by Nordlander et al. (2012). 
Here, we employ the photometric effective 
temperatures deduced from two different em- 
pirical colour calibrations of Casagrande et al. 
(2010) and Alonso et al. (1999) for dwarfs 
and those from calibrations of Alonso et al. 
(1996) for giants. Following Nordlander et al. 
(2012), we consider the SGB stars to be- 
long to the dwarf calibration. It is worth not- 
ing that the calibrations of Casagrande et al. 



(2010) are also valid for subgiants. The ac- 
curate determination of surface gravity of the 
GC stars takes advantage of their common dis- 
tance that can be evaluated from the compar- 
ison of the observed and theoretical colour- 
magnitude diagrams. The gravity was com- 
puted from the usual relation between log g 
and L, T e g, and mass of the observed stars. The 
used stellar parameters are indicated in Table 1 . 
Microturbulence velocities were derived from 
Fe n lines by Korn et al. (2007). 

The element abundances determined from 
individual iron lines turn out very similar for 
each pair of the SGB and RGB stars. As in 
HD 122563, the Fei lines arising from the 
ground term give higher abundances compared 
to those from the other Fe i lines, by 0.18 dex 
and 0.20 dex, on average, for the SGB and 
RGB stars, respectively. Therefore, mean Fe i 
based abundances were computed ignoring the 
low-excitation (E exc < 2 eV) lines of Fei. In 
contrast to HD 122563, we find differences of 
0.06 dex and 0.04 dex for the SGB and RGB 
stars, respectively, in mean abundance derived 
from the Fe n lines between applying gf- values 
from RUm and MB09. This is due to the 
smaller number of lines measured in the GC 
stars. Table 1 presents the obtained results. 

We find that the NLTE (S H = 0.1) abun- 
dances determined from two ionization stages 
are consistent within the error bars in the RGB 
stars, when using the RUm ^/-values, with 
(Fei - Fe ii)nlte = 0.07 dex and 0.02 dex 
for N11093 and N14592, respectively. In 
LTE, the corresponding differences amount to 
-0.1 1 dex and -0.16 dex. For the SGB stars, a 
discrepancy in NLTE abundances between Fe i 
and Fe n does not exceed 1 <r when employing 
the r e ff scale of Alonso et al. (1999) and the 
RUm ^/-values, but increases significantly for 
the hotter temperature T e ff = 5978 K based on 
the Casagrande et al. (2010) calibrations. 

When applying the 3D corrections pre- 
sented by Korn et al. (2007), we come to 
large discrepancy between Fei and Fen not 
only in LTE, but also in NLTE, namely, (Fei 
- Fe ii)nlte+3d is equal -0.2 dex, on average, 
for the two RGB stars, -0.2 dex and -0.1 dex 
for the SGB stars on the scale of Alonso 
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Table 1. Iron abundances loge and abundance differences A = Fei - Fen for the stars in 
NGC6397 



Star 


K 


logg 


6 
kms 


LTE 




NLTE 




Fei 


Fen 


Fei 


Fen 


A 




IRFM2010 1 














N5281 


5978 


3.73 


1.75 


5. 31 ±0.07 (15) 


5.34 5 ±0.07 (6) 


5.52±0.07 


5.34 5 ±o.07 


0.18 


N8298 


5978 


3.73 


1.75 


5.31±o.06(15) 


5.36 5 ±0.03 (6) 


5.52±0.06 


5.36 5 ±o.03 


0.16 




IRFM 1999 2 














N5281 


5816 


3.69 


1.75 


5.19±0.08 


5.25 4 ±o.09 


5.38±o.08 


5.25 4 ±o.09 


0.13 


N5281 


5816 


3.69 


1.75 




5.31 5 ±0.07 




5.31 5 ±0.07 


0.07 


N8298 


5816 


3.69 


1.75 


5.19±o.oe 


5.27 4 ±o.09 


5.39±o.06 


5.27 4 ±o.09 


0.12 


N8298 


5816 


3.69 


1.75 




5.33 5 ±o.03 




5.33 5 ±o.03 


0.06 




IRFM 1996 3 














N11093 


5100 


2.51 


1.7 


5.09±o.io (15) 


5. 16 4 ±0.10 (9) 


5.27±o.io 


5.16 4 ±o.io 


0.11 


N11093 


5100 


2.51 


1.7 




5.20 5 ±o.06 




5.20 5 ±o.06 


0.07 


N 14592 


5120 


2.58 


1.6 


5.09±o.08 (15) 


5.21 4 ±o.09 (9) 


5.27±o.09 


5.21 4 ±0.09 


0.06 


N 14592 


5120 


2.58 


1.6 




5.25 5 ±0.05 




5.25 5 ±o.05 


0.02 



Notes. Numbers in brackets indicate the total number of lines used in calculating the mean. 

1 r eff is based on calibrations of Casagrande et al. (2010) for dwarfs and subgiants, 

2 r ef f is based on calibrations of Alonso et al. (1999) for dwarfs, 

3 T tf! is based on calibrations of Alonso et al. (1996) for giants, 

4 based on gf- values of MB09; 5 based on g/-values of RUm. 



et al. (1999) and Casagrande et al. (2010), re- 
spectively. 

The above results may cast a shadow of 
doubt on 3D calculations. However, our inter- 
pretation is that one cannot simply add 1D- 
NLTE and 3D-LTE results when both NLTE 
and 3D effects are significant. Exactly this is 
the case in metal-poor stars. 

4. Centre-to-limb variation of solar 
Fe I lines 

At present, one cannot perform full 3D and 
NLTE line-formation calculations for iron due 
to its complex atomic structure. In this section, 
we investigate whether NLTE calculations cou- 
pled to average (3D) models can reproduce ob- 
servations of the centre-to-limb variation of the 
solar Fei 7780 A and 6151 A lines. The mea- 
sured equivalent widths (EW) were taken from 
Pereira et al. (2009). We apply the C0 5 BOLD 
solar 3D and average (3D) models and also 



the MAFAGS-OS solar ID model. When com- 
puting the theoretical EW, the iron abundance 
has been adjusted so that the models fit the 
observations at disk-centre. For the (3D) and 
ID model atmospheres, all the line profiles 
were computed with a microturbulence of = 
0.9 kms" 1 . 

The obtained results are presented in Fig. 3. 
It is evident that the classical model atmo- 
sphere fails to reproduce the observations. We 
do not show the NLTE results for the ID model 
because NLTE leads to weakened Fei lines 
and larger discrepancy between the theory and 
observations. The NLTE calculations with the 
(3D) model reproduce nearly perfectly the be- 
havior of Fei 7780 A, but the (3D) model is 
unsuccessful in the case of Fei 6151 A, both 
in LTE and NLTE. Full 3D-LTE line-formation 
calculations correctly predict the growth of 
£W(Fei 6151 A) toward the limb, albeit with 
overestimated theoretical EWs. The latter can 
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Fig. 3. Solar observed (rilled circles) and predicted centre-to-limb variations of line equivalent widths for 
Fe 1 7780 A (left panel) and Fe i 615 1 A (right panel). Filled rombs (only for Fe i 615 1 A) correspond to full 
3D-LTE line formation, crossed open circles to the NLTE calculations with the (3D) model, open circles to 
LTE with (3D), and open rombs to LTE with the MAFAGS-OS model. 



be due to neglecting the departures from LTE 
for Fe i. 

These results illustrate the varying sensitiv- 
ity of different Fe i lines to 3D effects. The Fe i 
7780 A line seems to be mostly affected by the 
change in mean atmospheric structure of 3D 
compared with ID model, while Fei 6151 A is 
more sensitive to atmospheric inhomogeneity. 

5. Test calculations with metal-poor 
average (3D) models 

It would be useful to produce a list of lines 
that are more sensitive to the change in mean 
atmospheric structure than to atmospheric in- 
homogeneity. Such lines could be reliably em- 
ployed in stellar abundance analyses with (3D) 
models. In this section, we check three lines of 
Fei and Fen 5018 A line in the full 3D, aver- 
age (3D), and reference ID models of VMP 
cool giant with T eB = 4480 K, log g = 2, 
and [M/H] = -3. The results are presented in 
Table 2 as the differences in abundance derived 
for a given line equivalent width between dif- 
ferent models. 

Despite the low iron abundance, all the se- 
lected lines are rather strong, with EW from 



70 mA to 95 mA. The full 3D-LTE line- 
formation calculations show that the effect of 
atmospheric inhomogeneity is weak only for 
Fe i 5 166 A (the row 3D - (3D) in Table 2, part 
I). In contrast, Fei 5232 A and Fen 5018 A 
are mostly affected by horizontal tempera- 
ture/pressure fluctuations. To examine a depen- 
dence of the 3D effects on the line strength, 
the calculations were also performed with gf- 
values reduced by one order of magnitude 
(Table 2, part II). The results differ dramati- 
cally: 3D effects are significantly reduced for 
Fei 5232 A and Fen 5018 A and, in contrast, 
are enhanced for the low-excitation lines, Fe i 
5 166 A and Fei 4994 A. In the case labelled 
"weak lines", all lines are predominantly af- 
fected by atmospheric inhomogeneities (the 
row 3D - (3D) in Table 2, part II). An excep- 
tion is Fe i 5232 A, with its minor (3D - (3D)) 
correction. 

We also performed LTE and NLTE calcula- 
tions with various kinds of averaged 3D mod- 
els. The full 3D model was decomposed into 
12 groups depending on the emergent inten- 
sity. The fraction of surface area with given in- 
tensity defines the weight of the corresponding 
group. By spatially and temporally averaging 
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Table 2. Abundance differences (dex) between different 4480/2/- 3 models 
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I. "strong lines" 


EW(mA): 


80 


75 


95 


70 


3D - IDlhd 


LTE 


-0.22 


-0.04 


+0.29 


+0.50 


3D - (3D) 




-0.02 


+0.08 


+0.35 


+0.43 


<3D> 12 - 1D L hd 


NLTE 


0.00 


+0.06 


+0.15 


+0.10 




LTE 


-0.22 


-0.14 


-0.02 


+0.10 




NLTE 


+o rn 


+n 04 


+n 05 


+0 0? 




LTE 


-0.05 


-0.04 


-0.01 


+0.02 


(3D) - IDlhd 


NLTE 


-0.07 


-0.01 


+0.10 


+0.09 




LTE 


-0.16 


-0.09 


-0.01 


+0.08 


II. "weak lines" 


EW(mA): 


25 


20 


35 


30 


3D - IDlhd 


LTE 


-0.39 


-0.23 


+0.02 


+0.20 


3D - (3D) 




-0.30 


-0.20 


-0.02 


+0.13 


<3D> 12 - IDlhd 


NLTE 


+0.04 


+0.08 


+0.13 


+0.10 




LTE 


-0.15 


-0.08 


+0.04 


+0.10 


<3D> 12 - (3D) 


NLTE 


+0.03 


+0.04 


+0.05 


+0.02 




LTE 


-0.09 


-0.06 


-0.02 


+0.02 


(3D) - IDlhd 


NLTE 


+0.01 


+0.04 


+0.09 


+0.08 




LTE 


-0.06 


-0.01 


+0.06 


+0.09 



each group we obtained the 12 average group 
models. The theoretical emergent line fluxes 
from all the average group models were inte- 
grated using the group weights. The obtained 
results are hereafter referred to as (3D)i2- 

In Fig. 4, we compare the Fei 5 166 A line 
fluxes from different models in the "weak 
line" case. In LTE, the line is stronger in the 
(3D)i2 model than in the (3D) and ID ones. 
In NLTE, the 3D effects are weakened and 
the (3D) model approaches the (3D) 12 one. In 
the "strong line" case, the difference ((3D)i2 
- 1D L hd) in LTE abundance is identical to 
(3D - 1D L hd) for this line. The (3D) [2 model 
seems to simulate the 3D model better than the 
(3D) model for the lines affected mostly by the 
change in mean atmospheric structure. 

As for the remaining lines under the LTE 
assumption, both (3D)n and (3D) models are 
inadequate for Fe 11 5018 A and for Fe 1 5232 A 
in the "strong lines" case. The 3D effects are, in 
part, reproduced by the (3D)n and (3D) mod- 
els for Fe 1 4994 A. 



We find that the NLTE effects for each Fe 1 
line are stronger in the (3D)n model than (3D) 
and ID ones. Only in the "strong line" case of 
Fei 5166 A, where the (3D)i2 model is work- 
ing, the 3D effects are fully compensated by 
the NLTE effects, and no difference in NLTE 
abundances was obtained between (3D) 12 and 
ID models. 

We caution against simple adding the 1D- 
NLTE and 3D-LTE corrections to obtain the 
correction to the abundance derived from clas- 
sical 1D-LTE analysis. For example, for the 
"strong" Fei 5 166 A line in the 4480/2/-3 
model, we calculated (3D - ID) = -0.22 dex 
under the LTE assumption. The NLTE correc- 
tion computed with the ID model amounts to 
Anlte = log snlte - log slte = 0. 16 dex. Their 
adding results in a combined correction of 
-0.06 dex, while we see no difference in NLTE 
abundances between (3D)i2 and ID models. 
When applying A NL te = 0.24 dex as computed 
with the (3D) model, the combined NLTE and 
3D correction is +0.02 dex. Therefore, the use 
of the average 3D model can be helpful, how- 
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Fig. 4. Theoretical NLTE (top panel) and LTE 
(bottom panel) fluxes (in ergcrrr 2 s _1 Hz -1 ) for Fei 
5166 A from the <3D)i 2 (dotted curve), <3D> (con- 
tinuous curve), and 1D LH d (dashed curve) models 
4480/2/- 3. All the computations were performed 
with log gf = -5.2 andf, = 1.8 kms~'. 

ever, we stress, only when the 3D effects are 
mainly caused by the change in mean atmo- 
spheric structure. 

6. Conclusions 

We summarize our findings as follows. 

1 . To within the error bars, consistent NLTE 
abundances were obtained from the two 
ionization stages, Fe i and Fe n, for the stars 
with well determined stellar parameters 
HD 122563 and red giants in NGC6397 
using classical plane-parallel model atmo- 
spheres. 

2. For HD 122563 and the stars in NGC 6397, 
the Fe i lines arising from the ground term 
reveal approximately 0.2 dex higher abun- 
dances compared to those from other Fe i 
lines. The obtained discrepancy between 



lines of different excitation lines cannot be 
explained by using the 3D corrections of 
Collet et al. (2007) and Hayek et al. (201 1). 
The 3D corrections predicted for the low- 
excitation Fe i lines are markedly too large. 

3. The use of the average (3D) models is justi- 
fied for particular Fe i lines with the 3D ef- 
fects mainly caused by the change in mean 
atmospheric structure of the 3D model 
compared to the ID one. Our NLTE cal- 
culations reproduce well the centre-to-limb 
variation of the solar Fe i 7780 A line, but 
they are unsuccessful for Fe i 6151 A. The 
metal-poor (3D) model was found to be ad- 
equate for the "strong" Fe 1 5 166 A (E exc 

! line, but inadequite for the lines sensi- 
tive to atmospheric inhomogeneities, i.e., 
Fe i 5232 A, Fe 1 4994 A and Fe n 5018 A. 

4. For Fe i, departures from LTE are expected 
to be stronger in 3D models than in ID 
ones. 

5. Simply adding 1D-NLTE and 3D-LTE 
corrections to the abundances determined 
from classical LTE analysis does not pro- 
vide 3D-NLTE results. 

Consistent 3D and NLTE line formation 
modelling is highly desirable for iron lines in 
metal-poor stars. 
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